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Abstract 

We consider a system of nonlinear equations that extends the 
Maxwell theory. It was pointed out in a previous paper that sym- 
metric solutions of these equations display properties characteristic 
of magnetic oscillations. In this paper I study a discrete model of 
the equations in two dimensions. This leads to the discovery of a 
new mechanism of vortex lattice formation. Namely, when a parame- 
ter corresponding to a magnetic field normal to the surface increases 
above a certain critical level, the trivial uniform-magnetic-field solu- 
tion becomes, in a certain sense, unstable and a periodic vortex lattice 
solution emerges. The discrete vortex solutions are proven to exist, 
and can also be found numerically with high accuracy. Description 
of magnetic vortices given by the equations is optical in spirit, and 
may be particularly attractive in the context of high-T c superconduc- 
tivity and the quantum Hall effects. Moreover, analysis of parameters 
involved in the discrete theory suggests existence of continuous do- 
main solutions — a conjecture that seems unobvious on grounds of the 
current topological and variational methods. 
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1 The proposed physical model 



Let A be the electromagnetic vector potential, and let the corresponding 
electromagnetic field be denoted by F A = dA. Further, let / be a real valued 
function. Consider the following system of equations 

dF A = (1) 



S(JF A ) = (2) 

- Af + \F A \ 2 f = vf. (3) 

The main goal of this paper is to describe a vortex-lattice type solution 
of this system of equations in two dimensions (cf. figure) by means of a 
finite difference approach. Results obtained in this way suggest that in the 
continuous domain limit such a solution consists of a continuous function / 
and a connection A whose curvature 2-form F A is also continuous. Moreover, 
the equations are satisfied in the classical sense almost everywhere. Note that 
the displayed numerical solutions are not defined at the vortex points. (Here, 
the horizontal axes are indexed by discrete lattice points.) Moreover, their 
scaling should be matched to a given physical system. 

As explained in Q (and in a simplified Abelian version in ||), the par- 
ticular form of this system is suggested by the geometry of principal U(l)- 
bundles. Purely phenomenologically dH)-© can be motivated by evoking 
heuristics used frequently in the optical literature. Namely, it is often as- 
sumed in optics that the nonlinearities arising in the interaction of radiation 
with matter can be accounted for by representing / as a series in tensorial 
powers of F A with coefficients characteristic of the material. Subsequently, 
one attempts to deduce properties of the coefficients from a microscopic the- 
ory. The shift of paradigm in this paper consists in assuming that / depends 
in a geometrically invariant way on F A via equation (H), and / remains es- 
sentially independent of the material except for a simple scaling. We will see 
below that in contrast to the soft optical nonlinearities, ([]])- (0) cannot be 
understood as a small perturbation of a linear system. However, interaction 
of radiation with matter is described from the "point of view" of the former, 
and the electronic processes inside matter are never discussed directly. 
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In order to explain superconductivity by means of a microscopic theory, 
one needs to display a mechanism that will let fermions overcome the obsta- 
cle imposed by the Pauli exclusion principle. A solution given by the famous 
BCS theory is based on the observation that when thermic noise is sufficiently 
low, it is energetically favored for electrons to join in pairs, known as Cooper 
pairs, which behave like bosons. The BCS theory is in a certain correspon- 
dence to the Ginzburg-Landau equations. These are nonlinear equations for 
a complex valued function if>, often interpreted as an order parameter govern- 
ing collective behavior of Cooper pairs. Periodic solutions of these equations 
in the form of vortices have been found by Abrikosov in |J. It must be em- 
phasized that strict validity of the BCS / GL approach, at least in its classical 
s-wave pairing version, is limited to low temperature superconductivity of 
metals. 

On the other hand, type II superconductivity is known to occur in ma- 
terials structurally different from metals, like YBCO, and at relatively high 
critical temperatures. As many researchers pointed out, this suggests that 
mechanisms other than those encompassed by the BCS theory may be re- 
sponsible for high temperature superconductivity. Those aspects of solid 
state theory that go beyond BCS seem particularly attractive in terms of the 
possibility of merging with the nonlinear Maxwell equations. It is possible 
that the new mathematical pattern introduced in this paper will be helpful in 
the description of the interaction of magnetic fields with composite particles. 
For illustration, consider the proposition that / describes a locally varying 
filling factor. In this interpretation, a part of the field gets entrapped in 
composite bosons, composite fermions, and Laughlin quasiparticles, which 
in turn see only the remaining y-fraction of the field. If / is a constant, this 
allows one to replace an electron picture with a composite particle picture, a 
suggestion that was present in science already ten years ago. However, if / 
is a vortex-type solution, the composite particles will feel the vortex in Fa, 
which should induce Josephson-type effects. Thus, while microscopic theory 
is always constructed with an a priori fixed filling factor, ([!])- (Q) would reflect 
the behavior of magnetic fields on a coarser scale. 

The above is meant to evoke some of the basic notions and ideas present 
in modern materials science. More thorough reviews can be found in the 



articles featured in the very incomplete list of references below (@-[12|) 
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2 Mathematics of a finite difference approxi- 
mation 



Consider the system ([!])- © on a two-dimensional flat torus T 2 . In this case 
S(fF) = implies *F = j for a constant B. In addition, dF = and 
F is the curvature of a certain connection A, provided its cohomology class 
satisfies [F] G 2ttZ. Thus, the system of equations reduces to 

"A/+y=I//, (4) 

and 

f B 

l — dV = 2ixK for an integer K. (5) 
M f 

Suppose / is a solution of the first equation with parameter B. Then for 
any c > 0, the function cf is a solution of (HI) with I? replaced by c5. At 
the same time, this rescaling does not affect (|5|) in any way, since the ratio 
B/f remains fixed. However, the system behaves differently with respect 
to rescaling of the independent variable. Indeed, suppose / satisfies both 
equations with parameters B, u, and N, then defining g(x,y) = f(cx,cy), 
we have that g satisfies (H) and (|5]) with parameters Be, z/c 2 , and N/c, while 
its period is 1/c in both directions. This is consistent with the experimental 
fact that as a magnetic field normal to the surface increases, vortices should 
eventually collide with one another. (In physical reality they will at that 
point disappear together with the superconducting state). This fact is also 
important mathematically, as it allows us to first obtain a solution of (Bj) with, 
say, f T 2 jdV < 2n, and then rescale the independent variable to satisfy 

It appears that the system (§])-© does not subdue itself to the standard 
techniques of variational calculus or topological analysis. In particular, per- 
turbative methods do not apply to (HJ), and no solutions arise as a result of 
bifurcation. In fact, in view of the theorem below and the results of numerical 
simulations, solutions are objects very unlike the familiar vortex solutions of 
nonlinear PDEs. To give some indication of the difficulties involved, consider 
the following. The equation (Q) is the Euler-Lagrange equation for the La- 
grangian L(f) — (| / |V/| 2 + B 2 /ln(/))/(/ f 2 ). However, this functional is 
neither bounded below nor above. Indeed, let us note that L(c n ) — > — oo for 
constants c n — > 0. On the other hand, let f n (x,y) = 1 + e + cos(2nnx). Since 
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haf n (x, y) > \ne, one easily checks that in this case L(f n ) — > oo. Therefore, 
the best we can expect is to discover local extrema. Additional difficulty 
stems from the fact that since ([|) always admits a trivial constant solution, 
we must devise a method of telling trivial and nontrivial solutions apart. 

We will now consider a finite difference model of the system (|)-(|). It 
is proven below that non-constant solutions of the discrete problem exist. 
The proof is independent of the number of points in the discretization (n 2 ) 
but relies on finite-dimensionality essentially, and does not admit a direct 
generalization to the analytic case. However, all the universal parameters 
used in the proof, like the L2-norm of / and B, are asymptoticly independent 
of n. Thus, we conjecture existence of the continuous domain solutions of (U) 
that satisfy the equation a.e. in the classical sense and retain the particular 
vortex morphology. 

It is convenient to introduce the following notation. / — I = 

n i j o 

~2 2 2 ; where indices i,j run through the discrete n-by-n lattice. Also, A 

denotes the common five-point periodic discrete Laplacian, i.e. A/(-, -) = 

^ (/(*£. J) + /& + J) + /& - 4/( J, £)) and V = (|, |) 

is the simplest two-point periodic gradient, i.e. say ^) = n(f(—, ^) — 

/(-, £)). In particular, the discrete integration-by-parts formula holds, i.e. 
— J(Af)g = /(V/, Vp). Consider the function 

$(/) = i/|V/| 2 + B 2 /ln(/). 

Pick arbitrary real numbers a, 6, c, fix a point (xq, t/q) = (— , — ), and a number 




Two submanifolds in R n , 

^, 6jC = {/ > : J f = a 2 , J - < i min/ = /(z ,2/o) = m n } (7) 
and its boundary 

dD n aAc = {/ > : j/ 2 = a 2 , 1 1 = ± min/ = /(x ,y ) = (8) 
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play a fundamental role in understanding the nature of critical points of 
$. Depending on the particular value of a, b, and c, the set D" 6c is either 
empty, an (n 2 — 2)-dimensional spherical disk, or it degenerates to a point. 
Consider the hyper-plane H n — {f : f(x Q , y ) = m„}. D™ b c is a spherical disk 
immersed in H n precisely when V, the point closest to the origin of an open 
submanifold given by f 4 = r — ^ — , is located inside the ball f f 2 < a 2 — ^f. 

o J o 

One easily finds V(x, y) = const = {n 2 — l)c for all x ^ x , y ^ y , and for it 
to be inside the ball, it is necessary and sufficient that a and c satisfy 

in2 - l) \ 2 <a 2 < (9) 



n 2 ri 2 



Conversely, if condition @ holds then D™ bc is a nonempty spherical disk. 
Formally, one needs to check in addition that m n is indeed a minimum of ev- 
ery function that satisfies the two other conditions in (^), but this is straight- 
forward. At this point I would like to point out that in the description of 
local minima of $ below, the parameter a is physical, and with good faith 
can be regarded as the L 2 -norm of the critical point, whereas the parameters 
b, c are auxiliary and will converge to as the density of discretization n 
increases to infinity. In particular, interpreting | as an approximate value of 
the integral of the reciprocal of the function where $ attains its local mini- 
mum is erroneous since the function develops a singularity at (xo,yo). The 
following theorem will be proven. 

Theorem 1 Fix constants B and a as above. For a certain choice of the 
constants b = b(n,a), and c = c(n,a), the function f — ► $(/) assumes local 
relative minima in D™ bc . In particular, the corresponding critical point, say 
fo > 0, satisfies the finite difference version of ($) everywhere except one 
point, i.e. 

B 2 

-Af (x,y) + — - = uf (x,y). for all (x, y) ^ (x , y ). (10) 

fo{x,y) 

Moreover, if B is sufficiently large, then fo is not a constant function. 

Proof. It is convenient to first treat D™ bc formally, and check that condition 
@ is satisfied later. Proceeding formally, let N denote an outward normal 
vector to dD™ bc inside its ambient sphere, i.e. iV is tangent to the (n 2 — 2)- 
dimensional sphere S a — {/ : J f 2 = a 2 , f(x ,y ) = m n } and points away 
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from the region D™ bc . The main task is to show that iV$ > 0. It will then 
follow from smoothness of $ (V is a linear operator and / — > ln(f) is smooth 
for / > 0) that it assumes a local minimum inside D™ bc . Let us choose for 
N the vector field defined by 

NAx v) = { ^ f ^ y) ~lik? for (^v)^(^,yo) (n) 
A ,yj \ otherwise, v ; 

where / G dD™ bc . A direct calculation shows that 

W) = i(-/ /A/+B2(1 -^) + /? A/ - B2 /?' 

It remains to analyze terms one by one. First, since a function in D™ bc 
assumes its minimum at (x , yo), we obtain 

-JfAf = J\Vf\ 2 + ^A f(x , y )m n > J | V/| 2 > 0. (12) 

o 

Next, by definition of D™ bc and (^) we obtain 

1 r 1 r 1 11 r 1 1 1 

6 J f J f n 2 m n J f b n 2 c 

so that F F 7 = an d therefore 

f(x,y)>c for 7^ (x ,y )- (13) 

As an immediate application we obtain 

u 

Finally, we obtain the following inequality 

o J 

= J- T T — ( f(i^ i) + f(i- i^) + f(— i) + f i^O - 4f fi 

> F F — i — f-4f(^ = -4 F F — 1 — 
= _4n 2 (i-^^) = 

(15) 
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Together, estimates (0), (|IJ), (|IJ) yield 

iW)>(l- (16, 

So far no assumption has been made about the constants. Now, in order 
to guarantee existence of local minima one needs to ensure that D™ bc is 
nonempty, and that the outer derivative is positive. Both these requirements 
are met if we pick 

a 

c = c(n) = — 2 , 

and 

b = bin) = -. 

With this choice of c, @ holds and if / satisfying the first two conditions in 
(|7|), assumed at some other point a value at least as small as m n , we would 
have f 4 > = \ K > \. which is a contradiction. Thus D™ hr is 

o J 

nonempty. On the other hand, inequality ([16]) becomes Nf&(f) > ^-(1 — 
-L^n 6 — -n 2 which implies Nf<&(f) > for n sufficiently large. Consequently, 
$ assumes a local minimum inside D™ bc . Equation (|T0| ) is automatically 
satisfied because it expresses the fact that the component of the derivative of 
$ which is tangent to the sphere / f 2 = a 2 vanishes in all directions except 
possibly the (xo, yo) -direction. 

Still, the local minimum / could a priori be a function constant every- 
where, except at the discontinuity in (x , yo)- This can be avoided by taking 
B sufficiently large. As B increases, the solution which is constant except at 
(x ,y ), say 

re : (x,y) ^ (x Q ,y Q ) 
m„ : otherwise 



fi{x,y) 



becomes unstable, i.e. it corresponds to a saddle point on the graph of $. 
Indeed, one checks directly that 

for <p tangent to the origin-centered sphere at fi, so that / f\<j) = 0. Now 
consider to be a non-constant eigenfunction of the discrete Laplacian on a 
torus and let A denote the corresponding eigenvalue. Shifting it if necessary, 
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we can assume that (p(xo,yo) = 0, so that <fi is orthogonal to f±. We have 
that / |V0| 2 = A/0 2 , and thus + ecj)) < only if § > A. Thus, /i 

is not a local minimum for B sufficiently large. Moreover, since / f 2 = a 2 , k 
converges to a as n increases. Naturally, one can make sure that A remains 
fixed regardless of n by always picking to be a discretization of the same 
trigonometric function. This shows that a choice of B which guarantees that 
fo is a nontrivial solution does not depend on the discretization. □ 

3 Closing remarks 

The proof above depends essentially on the fact that all functions and mani- 
folds are discrete. Indeed, the constants b = b(n), c = c(n) we have used tend 
to infinity as n — > oo, and some of the estimates make no sense in the limit. 
However, it is important that the magnetic induction B and the L 2 -norm of 
/ do not depend on n. On the other hand, simulation shows that the dis- 
crete solution /o retains its particular morphology independently of n, and is 
always subharmonic. In addition, multiplying (|10D by fX{(x,y)^(x ,vo)} an d in- 
specting the vicinity of the singularity we easily obtain v ~ [B 2 + J | V/| 2 ), 
which is also expected to converge. In summary, we have enough evidence to 
believe that the figure above shows a good approximation to a strong solution 
of the continuous version of equation (^) that would posess a vortex of Lip- 
schitz regularity. It is also consistent with the one-dimensional case, where 
solutions of the analogous nonlinear equation can be expressed in terms of a 
closed- form integral. 

I point out for completeness that continuity of such a positive solution f 
guarantees continuity of Fa = and by rescaling the independent variables 
we can satisfy the condition that the cycle of Fa be an integral multiple of 
2tt. This is sufficient to solve dA = Fa for A on a compact surface, and 
the solution A has sufficient regularity to retain its geometric interpretation 
as a connection 1-form. On the other hand, A contains all the information 
necessary to derive the basic tenets of superconductive electronics, like the 
Josephson effect. 

It is important to realize that none of the propositions stated above nec- 
essarily rely on the lattice being a simple square lattice. Most likely, a hexag- 
onal lattice setting would yield the same qualitative results. The only reason 
for using a square lattice is to avoid overwhelming numerical complexity in 
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experiments, as well as arithmetical nuisance in theory. 

I introduced the system (HI)-© i n 1993, and later investigated its prop- 
erties in my thesis. Some of those early results are contained in jTJ. The 
geometry was first given a loose but essentially correct physical interpreta- 
tion in 0. 

This work is a departure from the bulk of current research I conduct in 
collaboration with Professor R. R. Coifman in the area of Computational 
Harmonic Analysis. Working with Raphy is a stimulating adventure, and I 
want to thank him here for enhancing my understanding of mathematics. I 
also wish to thank my friend Fred Warner for proofreading the manuscript. 
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